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Several analyses of the microwave sky maps from the Wilkinson Microwave Anisotropy Probe 
(WMAP) have drawn attention to alignments amongst the low-order multipoles. Amongst the 
various possible explanations, an effect of cosmic topology has been invoked by several authors. We 
focus on an alignment of the first four multipoles {£ = 2 to 5) found by Land and Magueijo (2005), 
and investigate the distribution of their alignment statistic for a set of simulated cosmic microwave 
background maps for cosmologies with slab-like topology. We find that this topology does offer a 
modest increase in the probability of the observed value, but that even for the smallest topology 
considered the probability of the observed value remains below one percent. 
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I. INTRODUCTION 

Several recent analyses of the WMAP satellite maps 
have pointed out an unexpected degree of alignment be¬ 
tween the low-order multipoles of the cosmic microwave 
background (CMB) anisotropy Uillliiillill 
Various explanations have been put forward for these 
alignments, ranging from statistical fluke or foreground 
contamination through to a genuinely cosmological inter¬ 
pretation in terms of breakdown of statistical isotropy. 
Such a breakdown would be a natural consequence of the 
Universe possessing a non-trivial topology of character¬ 
istic scale comparable to the observable Universe (for a 
selection of cosmic topology review papers see Ref. M)- 

In this paper, we do not seek to address the interpre¬ 
tation of the observational data, but rather aim to test 
whether or not slab-space cosmic topologies give rise to 
the kind of alignments that are tentatively reported to 
have been observed in the first-year WMAP data. The 
observational indication is that there exists a preferred 
direction for the low multipoles. For instance, Tegmark 
et al. 01 and de Oliveira-Costa et al. 0 noted that the 
quadrupole and octupole were closely aligned with one 
another, and approximately planar. Land and Magueijo 
0 (hereafter LM) sought the alignment for each mul¬ 
tipole £ that maximized the proportion of power con¬ 
tributed by a single m mode, and noted that the align¬ 
ments of the first four multipoles were much closer than 
would be expected under statistical isotropy. These au¬ 
thors have all suggested that such alignments may be an 
indication of a slab topology where only one dimension 
is compact (finite and unbounded). 

The principal aim of this paper is to determine whether 
the LM alignment is a prediction of slab-space cosmic 
topology. We simulate CMB maps for spatially-flat slab 
topologies, for different sizes of the compact dimension, 
and derive the statistics of the alignments as defined by 
LM. We find that the degree of alignment in the observed 
data remain anomalous even in slab-space topologies. 


II. COSMIC TOPOLOGY AND THE CMB 

If the Universe has a non-trivial topology, that would 
lead to a breakdown of global isotropy. The spherical 
harmonic expansion coefficients aim of the observed map 
of CMB temperature anisotropies would then no longer 
be uncorrelated random variables, their correlation ma¬ 
trix having off-diagonal terms. Riazuelo et al. Pill 
found this correlation matrix for many multiconnected 
spaces, by computing the eigenmodes of the Laplacian 
with boundary conditions reflecting the particular topol¬ 
ogy 

Let us briefly recall how this is done. The power spec¬ 
trum coefficients, the Ce, are computed via the formula 

Ci (X J ejik)P{k)k^ dk, ( 1 ) 

where P{k) corresponds to the initial power spectrum of 
cosmological perturbations of wavenumber k, and Qi{k) 
are transfer functions defined in Ref. . 

In the case of Gaussian perturbations, the full statis¬ 
tical information is encoded within the two-point cor¬ 
relation function. The observable quantities correspond 
to the aim coefficients of the decomposition of the tem¬ 
perature field in spherical harmonics. The quantity 
we are interested in is therefore the correlation matrix 
{aimoiiim')- ^ simply-connected Universe, the cosmo¬ 
logical principle implies that the cosmological perturba¬ 
tions must be statistically isotropic. This in turn implies 
that the above correlation matrix is necessarily diagonal: 
{aimO-*iim') ~ where Ci is defined above. In 

the case of a multi-connected Universe, the Universe is no 
longer isotropic and the correlation matrix has non-zero 
off-diagonal components. 

In a spatially-flat multi-connected Universe, the eigen¬ 
modes of the Laplacian can be decomposed into the 
usual basis of spherical harmonics and spherical Bessel 
functions as 


T^s=£!;imk{kr)Yr{e,p) 


( 2 ) 
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where the index s distinguishes between modes with iden¬ 
tical wavenumber k. With this decomposition, the coef¬ 
ficients of the correlation matrix read 

cx ^ esimesV^*Qi(.k)Qi,{k)P{k). (3) 

k,s 

Once the coefficients are known, all the statistical 
information about a given topology can be computed us¬ 
ing an existing CMB code. This method therefore nat¬ 
urally takes into account all the contributions (Sachs- 
Wolfe, Doppler and Integrated Sachs-Wolfe (ISW)) to 
the temperature anisotropies, unlike the analytical for¬ 
mula of Ref. [2 which performs an estimate of the Sachs- 
Wolfe term only in the large-wavelength limit. 

Should one take into account the Sachs-Wolfe contri¬ 
bution only, then the CMB maps would exhibit sets of 
pairs of circles whose temperature pattern would per¬ 
fectly match. These correlated temperature patterns 
arise from the fact that we see two copies of the same re¬ 
gion along different lines of sight. The Doppler and ISW 
contributions reduce this correlation because the Doppler 
term depends on the direction in which the electron ve¬ 
locity field is observed, and the ISW effect depends on 
the photon history along the line of sight M- The simu¬ 
lated maps neglect reionization. As with the ISW effect, 
reionization tends to blur the topological signature in the 
correlation matrix, so that our simulated maps exhibit a 
stronger departure from statistical isotropy than more 
realistic maps. 

Following LM’s suggestion, we restrict our study to 
slab spaces. Our computational method requires all di¬ 
rections to be finite, so we chose rectangular tori with 
dimensions of the form 15 x 15 x X, labelled T[15,15, X] 
in the notation of Kunz et al. where X = 1,2,..., 15 
and the sizes are in Hubble radius units. Recall that the 
distance to the last-scattering region today is around 3.1 
Hubble radii in a flat ACDM model with ~ 0.7. The 
dimensions of size 15 Hubble units are essentially infinite, 
which we checked by comparing the correlation matrix of 
T[15,15,15] to that of a standard simply connected, in¬ 
finite universe. Hence T[15,15,A] is a computationally- 
favourable approximation to a slab space (that is, a space 
with only one compact dimension). Another reason to 
consider slab spaces is that matched circles searches have 
so far given negative results [l^, so that it seems likely 
that only topologies exhibiting a small number of circles, 
such as slab spaces, are compatible with the data. 


III. RESULTS 

A. The Land—Magueijo statistic and its observed 
value 

LM devised a statistic to study the alignments of mul¬ 
tipoles. For each multipole, they found the orientation of 
the coordinate axes which maximized the concentration 


of the multipole power into a single m, defining 


= max 


Clr, 


m.n (2£ + 1)C£ ’ 


( 4 ) 


where n is the coordinate axis orientation, Ci the usual 
power spectrum, and Cim measures the power at a single 
m value, defined as Ciq = ja^ol^ and = 2 la^ml^ for 
m > 0. The vector is defined as the one which provides 
the maximum value of the statistic ri for each £. 

They noticed a strong alignment of the of the lowest 
four multipoles, £ = 2,...,5 and quantified this by defining 
the mean angle between the six different pairs of align¬ 
ments. As the orientation vectors are headless (the same 
results are achieved by interchanging with —n^), one 
must choose the angle 9ij which is less than 90 degrees. 
The average alignment angle is then 


9 = mean (0^), i,j = 2,3,4, 5 with i ^ j ■ (5) 

They evaluated this alignment angle for the Tegmark, 
de Oliveira-Costa, Hamilton (TOH) cleaned and Wiener 
filtered maps of Ref. and found the values of 22.4° and 
22.3° respectively.^ We have confirmed this result using 
the code described below. 

For a Gaussian map, the predicted value of 9 is one 
radian. This arises as follows.^ Remembering that the 
alignment direction is a headless vector, the average sepa¬ 
ration between two such vectors can be found by rotating 
coordinates so that one is at the north pole. The other 
then has one end uniformly distributed in the northern 
hemisphere, so the average angle is given by the average 
distance of a point in the northern hemisphere from the 
north pole, which is one radian. 

Land and Magueijo found that the low observed value 
arose in only 5 out of 5000 simulated Gaussian maps, a 
result which we confirm below. On the face of it, this 
strongly excludes statistical isotropy, but one does need 
to bear in mind the strong a posteriori selection of their 
statistic; for instance the signal would be much weaker if 
the average angle included alignment with even just the 
sixth multipole. 

It is also found that because the statistic involves max¬ 
imization, it can be highly sensitive to small changes in 
the map, because multipole concentrations can have near 
double maxima for different m, giving completely differ¬ 
ent alignment directions. For example, while in the result 
above using the TOH maps, the quadrupole alignment 
was one which maximized the m = 2 multipole, with n 2 
in the direction {b,l) Ri (60, —100), we find (as did LM) 
that for the Lagrange internal linear combination map 
(LILC) produced by Eriksen et al. 0 a completely dif¬ 
ferent orientation is selected for the quadrupole, this time 


^ Their paper quotes ‘of order of 20 degrees’; we thank Kate Land 
for providing the actual value. 

^ Thanks to Kate Land for providing this argument, not given in 
their original paper. 
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maximizing the m = 0 multipole (in the TOH map this 
maximum is just slightly less than the one chosen), which 
significantly increases 6 to 55.2°. Similarly, in the case 
of the internal linear combination map (ILC) produced 
by the WMAP team da, while the preferred axis of the 
quadrupole is similar to that in the TOH maps, this time 
it is the direction of the octupole that is the cause of a 
discrepancy that increases 0. We note this as a caveat; 
in both these cases of near double maxima, three of four 
maps agree as regards the direction of the preferred axis 
for the concerned multipole. The maps differ in the de¬ 
tails of how foregrounds are removed from them. 

Again, we do not here seek to address data-related as¬ 
pects such as the likelihood of finding the observed 0 
given the uncertainties in foreground subtraction, the 
sky-cut, cosmic variance etc. While these issues are 
important to consider when determining whether the 
claimed detection of a breakdown of isotropy is signif¬ 
icant, papers that reported the detections in the first 
place have delved into such questions to some extent, 
and there is only limited progress one can make with an 
a posteriori detection. Aspects of fore grou nd contami¬ 
nation have been studied in Refs. [1.IRIlSL IT^ 1^ . For 
instance, Copi et al. Q considered the nature of known 
foregrounds without finding any clear connection to low-f 
alignments, at least as regards the LM kind. 

Despite the above-mentioned caveats, the value of 0 
found by LM is very low, and they speculate that it 
may be a signature of slab topology, already invoked in 
Refs. 0, 0 as a possible explanation of the quadrupole 
and octupole planarity. Our aim is to test this sugges¬ 
tion by evaluating the distribution of the 0 statistic for 
simulated slab topology maps, and checking whether, rel¬ 
ative to the case of a trivial topology, such topologies can 
better explain a low 0. 


B. The Land—Magueijo statistic for simulated 
maps of slab topology 

Given the full ensemble average correlation matrix 
{aimaiim') spherical harmonic coefficients for slab 

space topologies of form T[15,15,A'], we create corre¬ 
sponding random realizations of the atm- These are then 
rotated through a two-dimensional grid of galactic angles 
b and I (we used spacings of one degree in ea ch) . The aim 
transform under rotations as (see e.g. Ref. 

i 

aim= Y. (6) 

m'— 

where (a, /?, 7 ) are the Euler angles corresponding to the 
rotation, and 7 ) is part of the representation 

of the Wigner rotation matrix [^. Rotating the aim over 
the entire grid, the rotation (and associated m) that 
maximized Vi was recorded, for each i. Thus 9 was found 
for each simulation. This was repeated for an ensemble of 
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FIG. 1: The ensemble average alignment {9) and its standard 
error for the T[15,15, X] topologies, as a function of topology 
scale X. The horizontal line is the Ganssian prediction. 


universes so that the distribution of 0 was built up. The 
ensemble average value {0) was also found. This was then 
repeated for different dimensions X of the rectangular 
toroid T[15,15,A:]. 

Figure n shows the ensemble average values (d), to¬ 
gether with their standard error, against X. For Gaus¬ 
sian maps, the expected result of one radian is accurately 
recovered. We further see that the large-scale topology 
maps behave essentially as Gaussian maps. Only for 
X < 3 do we begin to see an effect of topology, with 
the ensemble average reducing slightly with respect to 
the Gaussian result. The shift is very small as compared 
to the observed value of 0. 

Although the shift of the ensemble mean is small, 
topology might nevertheless help explain the observation 
if it alters the distribution of 0 for small angles. Figure[21 
shows the full distribution of mean angle 0 obtained from 
10,000 realizations of the corresponding topologies. Such 
a large number of simulations was used in order to trace 
the tails of the distribution accurately. The value of one 
radian, which is the expected mean angle for a Gaussian 
random realization, is shown with a thick dashed line, 
and the observed value is shown as a dotted line. While 
the distribution of 0 extends a little towards smaller val¬ 
ues as the size of the smaller dimension of the toroid 
decreases, the observed value remains significantly low 
even for the smallest dimension considered. 

We quantify this further in Figure 0 which shows the 
fraction of the ensemble giving a value at least as low 
as the observed one. The uncertainties are estimated us¬ 
ing the Poisson error on the number of such ensemble 
members. We confirm the result of LM that approxi¬ 
mately 0.1% of Gaussian skies give the observed value or 
lower. For small topology scales we see an enhancement 
in this fraction, but even for the smallest topology con¬ 
sidered, the probability remains below one percent. We 
conclude that the observed alignment is not predicted by 





4 





Degrees 


FIG. 2: The full distribution of 0 for three T[15,15, X] topologies, X = 15,4,1, and for a Gaussian random realization. The 
dashed line shows the mean value and the dotted line the observed value. Note the extended low-end tail of the T[15,15,l] 
distribution, corresponding to the increased probability of this topology being able to produce the observed result. 
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FIG. 3: The fraction of realizations of a given T[15, 15, X] 
topology which provide a result at least as low as the observed 
result {9 = 22.5°), shown as a function of the characteristic 
topology scale X. 


slab-space topologies. 

We have only asked whether the topologies within our 
set can explain the result found by LM. This is not 
achieved even by the smallest topologies we consider. 
However we note further that such small topologies are 


almost certainly already excluded by other observations 
|22|. For instance, a harmonic space analysis of pre¬ 
cisely these same simulated topologies suggests that 
X < 3 is excluded by comparison to WMAP data (see 
also Ref. H^). Such topologies are also constrained by 
the null results of matched circles tests for such topolo¬ 
gies cmi, corresponding to a limit of about X ~ 5. 
[The very tentative indications of dodecahedral topology 
found in Ref. are interesting but of no significance 
for the topologies we consider here.] We have found that 
such topologies also do not receive any real support even 
from the observed alignments of the type discussed here. 


C. The Land—Magueijo statistic for ILC 
simulations 

Stepping somewhat aside from the main drive of this 
paper, we have made a small study of foreground effects 
by finding the distribution of 9 that results from the 
10000 ILC simulations provided by Eriksen et al. |^. 
These simulations additionally (over an assumed Gaus¬ 
sian CMB sky) contain the residual level of foregrounds 
that can be expected from the ILC method of foreground 
subtraction. The distribution of 6 that results is much 
broader and flatter towards smaller angles. That is ex- 
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pected because the simulations contain residuals of the 
galaxy, so that a preferred axis going through the galac¬ 
tic poles or thereabouts will be expected for the lower £ 
and hence 9 will be smaller. 

This result does not, however, explain the Land- 
Magueijo result, as the alignment they find (which is 
identified in the TOH maps) is not directed towards the 
galactic poles. It seems that if there were a preferred 
axis in a Gaussian CMB sky contaminated by galactic 
residuals, then it would be in the direction of the poles. 
This agrees with the conclusion of Copi et al. 0 that the 
alignments seen do not correlate with known foregrounds. 
Also, note that the real ILC maps do not show the same 
alignment direction effect, the effect being spoilt due to 
the previously-mentioned presence of a near double max¬ 
ima in one of the four multipoles considered. 

To summarize, foreground contamination could ex¬ 
plain a small 9 about the poles, but the observed ori¬ 
entation is not explained. Further, the alignment effect 
is detected in maps using a different foreground cleaning 
method P|. Hence this may be a case where a curious 
real feature is being obscured by the presence of galactic 
contamination (and cosmic variance). Interestingly, this 
point has already be made for related alignment statistics 
by Slosar and Seljak [^ . 


IV. CONCLUSIONS 

Our main results are as follows. We have confirmed the 
observed value of 9 found by Land and Magueijo, while 
noting that it is quite dependent on the choice of maps 
used. We have also confirmed their result that Gaussian 
skies have only about 0.1% chance of finding a value as 
low as that observed in the TOH maps. By analyzing a 
set of slab-topology maps, we have found that there is a 
slightly-enhanced probability of such a low value being 
obtained, but in absolute terms it remains extremely un¬ 
likely. We conclude that slab topology is not the expla¬ 
nation for the multipole alignment found by Land and 
Magueijo. The resolution must lie elsewhere, perhaps 
in other topologies, or instead in other cosmological as¬ 
sumptions, or in foreground or instrumental noise. 
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